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The magnetization process in the one-dimensional Kondo lattice model for the doped {uc < 
1) case is studied by the density matrix renormalization group method. A rapid increase of the 
magnetization is caused by the collapse of the intersite incommensurate correlation of / spins. On 
the contrary, the intrasite f-c singlet correlation survives in the larger magnetic field. The crossover 
from large to small Fermi surfaces for majority and minority spins is observed, whereas the Fermi 
surfaces are always contributed by / spins. A magnetization plateau appears with the magnitude 
of 1 — nc. Both ends of the plateau are related to the coherence temperature and the Kondo 
temperature which are characteristic energies essential in heavy electron systems. 
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In heavy electron systems, metamagnetism has been 
attracting great interest for a long time. A typical com- 
pound CeRu2Si2 shows metamagnetic behavior at the 
magnetic field Hm = 7.7 T In a de Haas- van Alphen 
(dHvA) experiment, the size of the Fermi surface changes 
at Hm H] . According to band structure calculations [g| , 
the measured Fermi surface is consistent with the itiner- 
ant /-electron band with the large Fermi surface below 
-ffivi- On the contrary, above Hm the measured Fermi 
surface is small, which is explained by the band struc- 
ture of LaRu2Si2; / electrons seem to be localized. In a 
neutron measurement, the incommensurate spin correla- 
tion disappears rapidly around Hm, but spin fluctuation 
on the single site survives above Hm 0- 

In the course of understanding above physics, it is a 
necessary step to exhibit a fundamental mechanism of 
the magnetization process in the basic model describing 
heavy electron systems. So far, numerical calculations 
for small clusters with a magnetization have been done 
in the periodic Anderson model (PAM) and the Kondo 
lattice model (KLM) |^-^. However, the unified picture 
of the magnetization process in the whole range of the 
magnetic field has not been established yet even in the 
basic models. 

The purpose of this Letter is to clarify the fundamen- 
tal mechanism of the magnetization process in the basic 
model on the basis of an accurate method: We calculate 
the magnetization, Fermi surfaces and correlation func- 
tions in the one-dimensional KLM with use of the density 
matrix renormalization group (DMRG) method |10|. We 
show that the magnetization plateau is related to two 
characteristic energies essential in heavy electron sys- 
tems. We also present the real-space picture and the 
band picture with respect to the quasi particle by com- 
paring the magnetization process of the KLM with that 
of the PAM. 

The one-dimensional KLM in a magnetic field is 



Sf 



(1) 



where cl^{cia) is a creation (annihilation) operator of a 
conduction electron on the i-th site {1 < i < L) with 
spin (T. S{ and are spin operators on the i-th site 
of the / spin and the conduction electron, respectively. 
The rescaled magnetic field denoted by h includes the 
(/-factor and the Bohr magneton fiB- The magnetization 
is defined by rn = nl — n| -I- — nj, where and 
are the densities of / spins and conduction electrons 
with spin a, respectively. The magnetization is related 
to S-fot = + 5f ) as TO = 2S^^JL and m takes 

the range of < to < 1 + Uc- 

We study the magnetization process of eq. (|l|) away 



from half filling (uc 



nj < 1). In the DMRG cal- 
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culation, we typically set the system size L = 40 under 
the open boundary condition and the number of states 
kept is taken up to 300. 

As for the doped case at /i = 0, a ferromagnetic metal- 
lic phase is known to exist in the large- J/i regime [ p^ . 
In this Letter we restrict to the paramagnetic metallic 
ground state without h. Figure |^ shows the magnetiza- 
tion process for t = 1 and J = 1 at ric = 4/5. A plateau 
appears with the magnitude of m = 1/5 as seen near the 
left end of Fig. |l|(a) (the enlargement will be given in 
Fig. H). The magnetic field of the left (right) end of the 
plateau is denoted by ho{hi). As ric approaches 1, hi 
increases and continuously connects to the spin gap Ag 
at half filling without h; Ag is associated with the char- 
acteristic energy scale called Kondo temperature Tk [0 . 
In order to illustrate the constituents of the magnetiza- 
tion, the average for the z-component of spin is shown in 
Fig. |l|(b). It turns out that the main component of the 
magnetization for < h < ho comes from / spins. The 
rapid increase of to beyond hi is still dominated by / 
spins and conduction electrons give the negative magne- 
tization; this is due to the internal field via the Kondo- 
coupling J and note that the density of states shifts to 
the direction opposite to that expected in this weak-h 
regime. On the contrary, the property for the higher 
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magnetic field up to the saturation /12 is characterized by 
conduction electrons. In Fig. |l|(c), the intra and intersite 
spin-spin correlations for the central i-th site are shown. 
At h — 0, the magnetic moments are mostly canceled 
by the incommensurate array of / spins, and the resid- 
ual moments are screened by formation of the collective 
Kondo singlet in the small- J/i regime [§[0. Even if h 
is switched on, the intra and intersite correlations shows 
minor change as long as < h < hi. However, once h 
exceeds hi, the intersite incommensurate correlation of / 
spins is destroyed abruptly to be the ferromagnetic one. 
On the contrary, the intrasite f-c singlet correlation sur- 
vives even in the larger magnetic field. Consequently, it 
turns out that the rapid increase of m beyond hi ^ Tk 
is caused by the collapse of the intersite incommensurate 
correlation of / spins. 




FIG. 1. (a) Magnetization curve and Fermi surfaces (b) 
Y.i{St)/L and Y,.{Sr)/L (c) intrasite correlation (Sf • S^) 
and intersite correlation (S,^ ■ S.^^j^), under the magnetic field 
h for t — 1 and J = 1 at Mc = 4/5. The inset in (a) is the 
enlargement of the small-Zi regime for kpi/iT and fepi/yr (see 
text). 

By the calculations for various Uc, we find that the 
magnetization plateau appears at m = 1 — ric . This can 
be intuitively understood if we take the real-space pic- 
ture which is relevant in the large- J/i limit: The number 
L(l — ric) is equivalent to the number of the hole whose 
site has an only / spin. For < h < hg, these / spins are 
polarized by the magnetic field, in order to avoid the en- 



ergy loss to destroy the intrasite f-c singlet correlation. 
As long as ho < h < hi, this state is still stable so that 
the magnetization plateau appears. Once h exceeds hi, 
the f-c singlet starts to be broken as shown above. 

We now turn to the location of the Fermi surface. The 
Fermi wave number kp^ for each spin defined as the sin- 
gular point in the momentum distribution function is de- 
tected by the period of the Friedel oscillation [|l3|. The 
open boundary condition used in the present calculation 
works as the perturbation which induces the density os- 
cillations: The one-particle charge density for conduction 
electrons under the open boundary condition is written 
as 



cos(2fcF|a;) . cos{2kFix) 



+ B 



cos[2(fcFt 



■ kYi)x\ 
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where (iV^(x)) = {N^{x)) + (iVf (x)) with {N^{x)) being 
the one-particle density with spin a The spin 

density (S"=^(x)) = ((A^f(a;)) - {Nl{x)))/2 has the same 
form as {N'^{x)), but with different coefficients A„ and 
B. 
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FIG. 2. One particle density in the real space (I) and its 
Fourier spectrum (II). (a) m = 0.2 and (b) m = 1.4 for t = 1 
and J = 2atnc=4/5. The Fourier transformation is carried 
out by using the central 30 sites. 



As shown in Fig. | 

and {SD = ((iV^) 



I, we calculate {Nf) = {N^^) -f (iVf^) 

^ - (iV^))/2 with {Nf,) = (etc.) 
for each magnetization. The Fermi wave number is de- 
rived from the peak structure in the Fourier spectrum 
of {Nf) and as in Fig. |ll. In Fig. |l(a) we see 

that the holes are located mutually keeping the interval 
of 4 sites in the case of m = 1/5 which corresponds to 
the plateau. This configuration is detected by the peak 
structure at 2(fcF| -t-fepj) = 27r/5 (mod 2it) in Fig. ^l(a), 
which persists in its position even if m changes. Addi- 
tionally, there appears the 2k-pa peak for each spin whose 
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position depends on m. As m increases (^1), the peak is 
detected more easily by (N^^^) rather than by {N^) and 
(S'f^). We see that the clear two peaks of 2kpi and ikpi 
appear and also that the amplitude for 2(fcF| + fcpi) be- 
comes quite small in Fig. ^l(b). We show the J = 2 
data in Figs. §(a) and @(b), since the J — 1 data con- 
tain the additional small peak which is not intrinsic as 
reported at ft, = |p^ . The resultant kp„ is shown in 
Fig. 0(a). At h — the Fermi wave number is given by 
kp^ = kpi = tt{1 + nc)/2 = 97r/10 [||: Namely, the 
Fermi surface is large, which includes the contribution 
from / spins. As h increases from to /iq, /cpt (^fj,) 
increases (decreases) to tt {iTnc = 47r/5), and they re- 
main the same as long as ho < h < hi. When h ex- 
ceeds hi, fcpT increases from and fcpx decreases from 
Trric. When h reaches the saturation field we obtain 
/cp| = TTUc — 47r/5 and fepj — 0, which are equal to the 
small Fermi surfaces of the completely polarized conduc- 
tion bands decoupled to / spins. 

We confirmed that the Fermi surfaces obtained with 
use of eq. (^ and the magnetization derived from the 
minimization of the total energy satisfy the following re- 
lation; fepcr/TT = (1 + Uc ± to)/2, where a =t (J,) cor- 
responds to +(— ) in the right hand side. This means 
that there is no phase transition and that the Luttinger's 
sum rule holds in a magnetic field. The expression of 
kpa is consistent with the argument on the basis of the 
generalized Lieb-Schultz-Mattis theorem The im- 

portant point is that the Fermi surface is always large, 
which involves the contribution from / spins, whereas 
the crossover from large to small Fermi surfaces is ob- 
served fl7|] . 

In order to analyze the crossover from large to small 
Fermi surfaces in detail, let us consider the PAM: 

HpAM (^Lc^+i^T + H.c.) +SfJ2 i^ia 

i(T ia 

+ ^E(/-^- + H-c.)+c/E"/T"a- (3) 

ia i 

The KLM is the effective model derived from the PAM 
in the regime of C/ ^ ■np{Q)V'^ under the symmetric con- 
dition £/ -I- U /2 = 0, where p{0) is the density of states 
of conduction electrons at the Fermi energy. Finite U 
makes the effective / level shifted up to the Fermi energy 
where the Kondo resonance responsible for the heavy 
quasi-particle band arises. From the viewpoint of the 
adiabatic continuation, we consider the simplest case of 
£/ = —U/2 = in eq. (||). In Fig. ||, the solid curve rep- 
resents the magnetization of the PAM with the Zeeman 
teYm.-hY,i{S(^+Sf) for the doped (n = nf+n, = 9/5) 
case. 

For h^ < h < h\ the plateau at m — 2 — n ap- 
pears as a consequence of the hybridization gap: A — 
-2t + 2y/WTV^. Here the critical fields h* {i = 0, 1) are 
given by 

h* = t {(-1)' - cos(n^) - i-iy^TTWW 



+ Vcos2(n7r) + (FA)2j, 



and the saturation field by hl+t[l+cos{mT)]. Namely, 
the plateau appears as long as the bottom of the up- 
per empty band with majority spin becomes lower than 
the top of the lower filled band with minority spin as 
h increases. In order to separate the contribution from 
conduction electrons, we also show the magnetization for 
V — 0, setting Uf = 1 and Uc — 4/5. Additionally, the 
discrete data of J = 1 and 2 in the KLM at ric = 4/5 
are shown, whose horizontal axis is scaled by hi. To fa- 
cilitate the comparison between the PAM and the KLM, 
the magnetization curve of the PAM is also scaled by K 
and we choose the hybridization V fixing t — 1 in eq. (|^) 
so that hi is equal to hi of the KLM. The V = case is 
scaled by the same unit of the corresponding V ^ case. 
Note that the data for the PAM are just located on those 
for the KLM for < ft < /ii. In Fig. |, the solid arrow 
indicates the magnetic field on which the sign of (S{ • S'i) 
changes in the KLM (see Fig.|l|(c)). 
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FIG. 3. Magnetization of the KLM for = 4/5 and of 
the PAM for n — 9/5, with t = 1. The dashed line represents 
the V — case of the PAM with Uf — 1 and ric — 4/5. Each 
data for V ^ and V = of the PAM is scaled by hi , where 
hi is set to be equal to hi of the KLM. The solid arrows show 
the point where {S{ ■ SI) of the KLM changes the sign. The 
open arrow indicates the crossover point from large to small 
Fermi surfaces (see text). 

As for the PAM, we see in Fig. ||that the magnetization 
for finite V approaches that for V — with Uf — 1 and 
ric = 4/5 as ft increases. Namely, the crossover from / to 
conduction electrons which gives the dominant contribu- 
tion to the magnetization as well as the Fermi surfaces 
occurs beyond ft*. Note that, even in the large-ft regime 
where the Fermi surface is dominated by conduction elec- 
trons, / electrons are coupled to conduction electrons; / 
electrons are always itinerant, but not localized. Since A 
is reduced to be the order of Tk = 2texp{-U/{8p{0)V^)} 
by finite U, the range where / electrons give the main 
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contribution to the magnetization process will be drasti- 
cally reduced by the many-body effects. As seen in the 
V ^ limit, the Fermi surface can be close to that of con- 
duction electrons beyond slightly larger h than hi ~ Tk 
as indicated by the open arrow. The difference of the size 
of the Fermi surface between V ^ and V = cases 
around the open arrow becomes small as n approaches 2 
(half filling) . Similarly, we expect that the crossover will 
be observed around hi ~ Tk in the KLM in the realis- 
tic J/t case near half filling, where the intrasite f-c 
singlet correlation remains. 

As seen in Fig. ^, the general behavior of the PAM 
is similar to that of the KLM. This implies that the 
fc-space (band) picture for the PAM continuously con- 
nects with the real-space picture for the KLM explained 
in Fig. 1^. Namely, for the PAM the smaller end of the 
plateau (/iq) is associated with the slope (i.e., velocity) 
of the lower hybridized band whose component is dom- 
inated by / electrons, and the larger one (hi) connects 
to the hybridization gap A as rt approaches to 2. The 
correspondence between the above picture and the result 
of the KLM leads us to find the following meaning of 
the magnetization plateau in the KLM: The smaller end 
of the plateau (hg) is associated with the velocity of the 
effective carriers whose main component is unscreened 
/ spins with density 1 — rij,, and the larger one (hi) is 
related to the intrasite f-c singlet correlation, which con- 
nects to the spin gap Ag as ric approaches 1. These two 
characteristic energies are considered to be related to the 
coherence temperature T* and the Kondo temperature 
Tk; the former and the latter are often observed as the 
characteristic temperatures which give T^ and logT be- 
haviors in the resistivity, respectively in heavy electron 
compounds. The quasi-particle picture with respect to 
the KLM obtained here via the response to the magnetic 
field is consistent with recent calculations of dynamical 
and finite-temperature quantities without h by N. Shi- 
bata, et al. jl^ 

When realistic factors which should enter are taken 
into account, the plateau tends to be smeared; the dif- 
ference of g factors between / and conduction electrons, 
anisotropy of the hybridization, the dispersion of the / 
band, etc. [|l9| Hence, the shape of m itself for < h < hi 
changes according to these factors. However, physics 
shown here offers significant pictures to understand the 
magnetization process of heavy electron compounds. As 
for CeRu2Si2, it was confirmed that there is no first order 
phase transition at Hm and the electronic state changes 
continuously This suggests that the crossover from 
large to small Fermi surfaces occurs with the /-electron 
number involved. Simultaneously, we pointed out that 
the localization of / electrons above Hm is not necessary 
to interpret both the results of the dHvA experiment |^] 
and the band structure calculations Furthermore, we 
show the rapid increase of the magnetization due to the 
collapse of the intersite correlation of / spins and the 
toughness of the intrasite f-c singlet correlation. These 



tendencies are observed by the neutron measurement . 

In conclusion, we show the fundamental mechanism of 
the magnetization process in the one-dimensional KLM 
and shed light on the heavy quasi-particle picture via the 
response to the magnetic field. 
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